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Abstract 

Recent studies suggest that the process of symmetry breaking after inflation typically 
occurs very fast, within a single oscillation of the symmetry-breaking field, due to the 
spinodal growth of its long-wave modes, otherwise known as 'tachyonic preheating'. We 
show how this sudden transition from the false to the true vacuum can induce a significant 
production of particles, bosons and fermions, coupled to the symmetry-breaking field. 
We find that this new mechanism of particle production in the early Universe may have 
interesting consequences for the origin of supermassive dark matter and the generation of 
the observed baryon asymmetry through leptogenesis. 
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1 Introduction 



Spontaneous symmetry breaking (SSB) is one of the basic ingredients of modern theories of 
elementary particles. It is usually assumed that SSB in Grand Unified and Electroweak theories 
took place in the early Universe through a thermal phase transition. However, it is also possible 
that some of these symmetries were broken at the end of a period of inflation |l| , when the Uni- 
verse had zero temperature and the negative mass term for the Higgs field appeared suddenly, 
i.e. in a time scale much shorter than the time required for SSB to occur. In this case, as was 
recently shown in Refs. 0, the process of symmetry breaking is extremely fast. The exponential 
growth of the Higgs quantum fluctuations is so efficient that SSB is typically completed within 
a single oscillation, while the field rolls down towards the minimum of its effective potential. 
This process, known as tachyonic preheating, leads to an almost instant conversion of the initial 
vacuum energy into classical waves of the scalar fields, in contrast with the process of 'para- 
metric preheating', in which the infiaton field performs many oscillations before reheating the 
Universe 0. 

In this letter we describe how this sudden transition from the false to the true vacuum 
can induce the non-adiabatic production of particles coupled to the Higgs. We also studied the 
consequences that this new process may have on the generation of the dark matter and the baryon 
asymmetry via leptogenesis. The phenomenon of particle production from symmetry breaking 
is analogous to the Schwinger mechanism where the role of the external electric field pulse 
is played here by the time-dependent expectation value of the Higgs field. It is also similar to 
the well known process of particle production by a time- dependent gravitational background 0, 
responsible for the observed anisotropics of the microwave background, as well as for Hawking 
radiation The difference with respect to the standard mechanism of gravitational particle 
production in an expanding Universe is that, in our case, quantum excitations of fields coupled to 
the Higgs are produced due to the sharp change in their masses during the process of symmetry 
breaking, instead of due to the quick growth of the scale factor. 

We will consider here a simplified model of SSB in which the Higgs instantly acquires a 
negative mass-squared term 0. This 'quench' approximation corresponds to the limiting case 
of a hybrid infiation model [[7[ satisfying the so-called 'waterfall' condition. We therefore assume 
that the complex symmetry breaking field starts in the false vacuum at the top of its potential 
V{(f)) = A(|0p — f^)^/4, with zero mean, (0) = 0, and initial conditions given by vacuum 
quantum fiuctuations. We will study the production of other fields, scalars x fermions ip, 
coupled to the Higgs with the usual scalar g'^\(f)\'^x'^ ^ind Yukawa h(l)il>%l) interactions. As we 
will show later, the backreaction of the produced x ^-nd ip modes on the evolution of the Higgs 
expectation value (vev) is negligible for the small g and h couplings we are considering, so that 
we can solve first the process of symmetry breaking, and then take the resulting evolution of 
the Higgs as a background field that induces particle production. 
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2 Symmetry Breaking 



The dynamics of symmetry breaking has been studied in detail in Refs. p|. Here we only 
summarise the main results needed for our analysis. At the initial stages of SSB, when re- 
scattering effects are still unimportant, the Higgs modes follow the linear equation 0^ + (/c^ — 
0fc = 0, where = Af^. With vacuum initial conditions, 0a:(O) = l/v^2fc, 0^(0) = 
— i/c0fc(O), all long- wavelength modes within the horizon {H < k < m) grow exponentially, 
0fc(^) = 0fc(O) exp(t\/ m? — kP'), while modes with k > m oscillate with constant amplitude. The 
exponential growth continues until the long-wave modes reach a value for which the effective 
Higgs mass becomes positive, i.e. when > v'^/V3, and the symmetry is broken soon after. 

We have chosen + | = \(j)l{t)(f)k{t)\ as a proper definition for the occupation numbers of 
the Higgs tachyonic modes [§]. This expression does not require an a priori definition of a mode 
frequency uJk, and matches smoothly the one used in for positive frequencies. For the growing 
modes, we have occupation numbers 
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that become exponentially large very quickly for the long wavelength modes, and drop abruptly 
for k > k^ = m(2mt)~^/^, which gives a natural ultraviolet cutoff for the problem. Note that, 
while the dynamics conserves (0) = 0, the Higgs dispersion also grows exponentially 

, m2 rmde m2 e^"^*(2mt-l) + l 

) = L ZZ^^ = 77^. 77^. e , (2) 



where this expression has been regularised, as described below. The time it takes for the system 
to break the symmetry, i.e. when (|0p(t*)) — f^, can be estimated as 

1, /327r% 

mt=^ ~ -log (— — 1 . (3) 



2 "V A 

This time depends only logarithmically on the Higgs self-coupling constant A. For typical values, 
A = 10~^ and v = 10~^Mp, the symmetry is broken within mt^, ~ 6, and the typical cut-off 
frequency becomes A;* ~ m/3. This means that, by that time, the occupation numbers of modes 
with k < k^, is exponentially large, 

ni^lQ ^ -e^"^** ~ ^ ~ 2 X 10^ (4) 
2 A 

These large occupation numbers allow us to treat these modes as semiclassical waves and match 
the solutions of the linear equations with the fully non-linear numerical lattice simulations 0]. 
The non-linear dynamics is studied by solving the real time evolution equations of classical 
fields, using a modified version of the lattice simulation program LATTICEEASY of Felder and 
Tkachev |]TU[. We start with initial fluctuations described by a Gaussian random field with zero 



mean, (0^) = 0, and regularised dispersion 



{mt)\u ^ mm - — . (s) 



2 



This prescription amounts to substituting quantum averages by ensemble averages and ensures 
that the physical masses and energies are not ultraviolet divergent. In Fig. |l], we show the result 
of the full non-linear evolution of the Higgs vev, and compare it with the approximate expression 

m ^ {mt)y/^ = |(l + tanh^^^^^) , (6) 

which ignores the strongly damped oscillations after symmetry breaking [Q. We have checked 
that these low-amplitude oscillations do not contribute to the non-adiabatic production of par- 
ticles; i.e. parametric preheating is inefficient after symmetry breaking, a result anticipated in 



Ref. [|T^ for the case of hybrid inflation. Note also that we are using the Higgs vev averaged over 
the lattice as a homogeneous background field, while it is actually a sum of tachyonic modes 
with different frequencies. However, in analogy with the generation of anisotropics during infia- 
tion, what drives the growth of long-wave modes k < k^, oi particles coupled to the Higgs is the 
coarse-grained Higgs vev over scales ~ which is practically the same as the average over 
the whole Hubble volume, thanks to the fact that all Higgs modes with H^^ < k < k^ grow 
essentially at the same speed. 



3 Particle Production 



We can now calculate the production of bosons and fermions coupled to the Higgs using the 
formalism of quantum fields in strong backgrounds [|l^, [l^. The mode equations in terms of 
rescaled fields, Xk{t) = a^^'^Xk and = a^^'^ip, 



are 



d^Xk+{k^ + ml{t)a^{t))Xk = 0, 
(i7^9^ -mF(t)a(t))^ = 0^ 



(7) 



where both the mass m{t) and the scale factor a{t) depend on time. In hybrid infiation models 
for which the quench approximation is valid, the rate of expansion is typically much smaller than 
the masses involved, and we can take the scale factor to be constant (a = 1) during SSB. We 
will only consider here the non-adiabatic production of particles due to the change of vacuum 



as it induces a sudden change in the inertia (masses) of bosons m| (t) 
m-p{t) = through the Higgs mechanism. 



(7^(101^) and fermions 



We have solved the mode equations (|^ and (H) both numerically, using the lattice Higgs 
vev evolution, and within the approximation (^), for which one has analytical solutions in terms 
of Hypergeometric functions [|12| . The number density of created particles, as seen by a future 



observer in the true vacuum, is then given by 



nx 



9s 



27r2a3 



dkk"^ \f3k\ 



(9) 



where X denotes either bosons or fermions; (3k are the Bogoliubov coefficients that relate the 
in {t — oo) and out {t +oo) mode functions; k is the comoving momentum, and we have 
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summed over spin indices {qs = 1 for scalars, 2 for spinors). In the case of charged fields, nx 
gives the number of particles, equal to that of antiparticles (i.e. gs = 2 for complex scalars, 2 
for Majorana fermions, 4 for Dirac fermions). 

Let us consider first the production of bosons. The mode functions Xk{t) of the scalar field 
are solutions of the oscillator equation (0) with time- dependent frequency ujI = k'^ + m^, and 
initial vacuum conditions, Xfc(O) = l/v^2cJfc and X[(0) = —iukXk. In this case, the Bogohubov 
coefficient at t ^ oo can be written as raf = |/5fcp = +u;^|Xfcp — ujk)/2uji^. The exact 

solution of the bosonic equation, 

XI + (^l{k) + ^(1 + tanh xf^ = , (10) 

with X = m{t — t*)/2, can be written in terms of Hypergeometric functions as [|l^ 

Xk = Ci{k) e'''\l + e™*)" F(a, 6, c; -e™*) + C2{k) e-*^*(l + e'"*)" F{a\ b*,c*; -e'"*) , (11) 

where r is a solution of the equation r(r — 1) + = 0, and 

a = T + i{uj^ — uj^) , h = T + + u^) , c = 1 + 2i ci;_ , (12) 

where uj± are the in/out asymptotic frequencies, uj-{k) = k and LO^{k) = yjk"^ + a^m?, and 
a = g/VX. These solutions match the initial conditions at the end of infiation, and can be 
analytically continued to t — > oo, in order to compute the boson occupation number = |/9fcp, 

B cosh[7rV4a2 ~ 1] + cosh[27r(co'+ — LU^)/m] 

sinh[27rcj_/m] sinh[27rcj+/m] ' ^ ' 

A similar analysis can be done in the case of fermions, where the first-order Dirac equation (H) 
for the two spinor components can be written as an oscillator equation with complex frequency. 



X'l + (fc2 + ml- im'p) Xk = 0, see Ref. p 



X'l + \^_{k) + —(1 + tanh xy - —(1 + tanh^ x)j = . (14) 

Given the initial vacuum conditions, Xfe(O) = (1 — mp j ujk)^!'^ and X^(0) = —icu^Xk, we can 
write the solution as in Eq. (pUf ), but with r = —ia. The fermionic occupation number can then 
be written as = = {uk — mp — ImXfcX^*)/2 cijfc. Substituting the t — > cxd solutions of 
the fermionic equations we find 

p cosh[27ra] — cosh[27r(ct;+ — c<j_)/m] 

Tlkia) = — ; — ; , (15) 

2 smh[27ru;_/mj smh[27ru;+/mj 

with u!- and uj^ the same as for bosons but with a = h/^/X. 

The occupation numbers obtained numerically from Eqs. (|^ and (|^), using the full non- 
linear lattice solution for the Higgs vev, agree very well with the analytical formulae (|13D and 
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([TB|), see Fig. 2, except at large momenta, where there is a small contribution from resonant 
particle production due to the strongly damped Higgs oscillations. 

We will now calculate the ratio of energy densities of the particles created at the end of 
symmetry breaking to the initial false vacuum energy density 

Po J 

where k = k/m and po = itl'^/4:X. We can obtain a fit to the final energy density of bosons and 
fermions as 

Pb/po - 2xlO-3(7,A/(a,1.3), (17) 
Pf/po - 1.5 X 10-3 (7,A/(«, 0.8), (18) 

where the fitting function is /(a, 7) = + 7^ — 7, and the coefficient a = g/\f\{h/\f\) for 
bosons (fermions), has been fitted in the range a G [0.01, 10]. 

We can also write the number densities for bosons and fermions, 

nB(a) ^ 1 X 10-3^,mV(a, 1.3)/a, (19) 
nF(a) ^ 3.6 X 10-^^,mV(",0.8)/a. (20) 

From these expression we see that the production of bosons and fermions from symmetry break- 
ing is more efficient when the Higgs mass and thus A is large, driving a very steep growth of 
the vev towards the true vacuum. In the case of very large couplings g"^, h? ^ A, the pro- 
duced particles are non-relativistic, Uk ~ rnx, and their energy density is given by px = rnyjiy^. 
Although the energy density increases with the mass of the particles produced, their number 
density saturates for large a to a value proportional to m^, and thus the number of particles 
produced cannot be made arbitrarily large by increasing g and h couplings. Furthermore, the 
physical momenta redshift as /cphys = k/a with the expansion of the Universe, and thus the 
number density of both relativistic and non-relativistic bosons or fermions decays as nx ~ a"^- 

Note that, unless the couplings are unnaturally large, the fractional energy density in bosons 
and fermions is always small, so we do not expect an important backreaction on the evolution 



of the Higgs condensate as the symmetry is broken. Moreover, contrary to the case of Ref. ||14|| , 
in which particles are produced long after symmetry breaking from non-linear rescattering, our 
mechanism of particle production from symmetry breaking gives an upper limit to the occupation 
numbers of bosons produced in the range H < k < m, even for arbitrarily large coupling g, 
which is of order < 1. This prevents us from using LATTICEEASY to compute their energy 
density and backreaction. 

We would like to explore now the cosmological consequences that this production of particles 
may have for the evolution of the Universe. 
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4 Leptogenesis 



Leptogenesis from preheating was first proposed in Ref. ||T5[ in the context of chaotic inflation, 
where the right-handed (RH) neutrinos were produced via parametric resonance through their 
couphng to the inflaton.Q Here we propose a novel scenario in which a large population of 
out-of-equilibrium RH neutrinos is produced in the process of symmetry breaking after hybrid 
inflation, via their coupling to the SSB field. 

Let us suppose that the symmetry-breaking field responsible for the end of hybrid inflation, 
and undergoing tachyonic preheating, is in fact the Higgs field of a grand unified theory (GUT) 
with (i? — L)-violating interactions, and that the same Higgs produces a significant population of 
right-handed Majorana neutrinos via the non-perturbative mechanism described above. These 
massive neutrinos decay (out of equilibrium) into leptons and standard model (SM) Higgses, 
generating a leptonic asymmetry that later gets converted into the baryon asymmetry of 
the Universe via sphaleron transitions JlF 



The Lagrangian terms relevant for leptogenesis are 

C=~hNI,<PNR-KNHlL, (21) 

where the first term determines the masses of RH neutrinos through the GUT Higgs mechanism, 
= h {(f)) = hv, while the second term determines the lepton- violating decays of the RH 
neutrinos into SM Higgses H and left-handed leptons II- Although both Yukawas h and hjy 
are actually 3x3 matrices, we will not consider here the complications of 3-generation neutrino 



masses and mixings, but will assume as usual [19| that only the lightest RH neutrino A^^i, with 



mass Ml, is responsible for leptogenesis. For analogous reasons, the CP- violating asymmetry e 
in the decay of the RH neutrino, which also depends on the structure of neutrino masses and 
mixing, will be taken free parameter. 

Assuming that the masses of the SM light neutrinos are generated by the see-saw mecha- 
nism m^. = {m £1171^2-,) a /Mi, the decay rate of the lightest RH neutrino can be written in 
terms of m^-^ and Mi as 

T.,-^^^n..M?...12.,0'k^ GeV, (22) 



Under these assumptions, the model has 5 parameters: the vev of the GUT Higgs v, the Higgs 
self-coupling A, the Yukawa coupling of the RH neutrino h, the mass of the light neutrino 
and the final reheating temperature Tj-h. We will now study the constraints that will determine 
the allowed range of parameters for a successful model of leptogenesis. 

The first requirement is that the RH neutrinos produced at symmetry breaking decay out of 
equilibrium. To ensure this we impose two conditions: i) that the Higgses produced via tachyonic 
preheating do not reach thermal equilibrium before the RH neutrinos decay; and ii) that the 

^Fermion production during preheating via parametric resonance was also studied in Ref. [ p^ . 
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Higgs-neutrino interaction rate is much smaller than the rate of expansion of the Universe, so 
that even if the Higgses thermalise, the RH neutrinos remain out of equilibrium. 

In order to satisfy condition i), we require that the Higgs self-interaction rate at the time 
of symmetry breaking is smaller than the rate of expansion, Fgeif < -f^sB = (27r/3)^/^ mf /Mp, 
where 

f (Pk ^ k 2\m 

^-'' = J (^^scif%(t)-=^j^^32^,^^^, (23) 

with (Tseif = and nl{t^) = g-fc^/afc^^ Note that, once ensured at symmetry breaking, this 
condition is always satisfied since nfp{t) ~ while H ~ t~^. This condition constrains the 
Higgs self-coupling to be 

ln(327r7A) V 6 ^ Mp ' ^ ^ 

Once chosen a GUT scale v, this constraint gives an upper bound on A, and thus on the maximum 
number of RH neutrinos produced in this model. We will take v = 10^^ GeV, and A = 10"'^, 
which satisfies the bound. 

The second constraint ii) requires Tint < H^b, where 



int 



d^fc ^ p k h?"m 

J7^ (Tint <(1 -nk) — = 21 2/00 2 /\^ ' (^5) 

[Ztt)'^ m 77"^ In (327r"^/A) 



with the Higgs- A?"! cross-section given by (Jint = iQ^2j^.ii , and where we have neglected the fermion 
contribution -C 1. This gives 

h ^ /27r\i/4 ^ V ^l/2 
M32^<i— ) [-^) ■ (26) 

For the chosen parameters, we can satisfy this constraint as long as h < 0.5. We will choose 
h = 0.01, i.e. Ml = 10^^ GeV. With these values of the parameters, the number density 
(|20|) and fractional energy density (|18|) of RH neutrinos at symmetry breaking becomes — 
1.4 X 10-^ and PnJpo ^ 7.4 x 10"^. 

We will now study the constraints associated with the time of decay of the RH neutrinos, 
ti = r^|. Since v, A and h have been determined by the previous constraints, these new 
conditions will give an allowed range for rriu-^^ . First of all, the Ni lifetime should be greater than 
the time of symmetry breaking, ti > t*, with t^, given in (^). For our parameters, we should 
satisfy 

We also have to ensure that the annihilation rate of A^^i into Higgses, Fann = CTann, with anni- 
hilation cross-section o"ann = /i'^/lGTT^Mf , is smaller than the decay rate, Fann < TtVi, otherwise 
no RH neutrinos would be left to produce the lepton asymmetry. This imposes the constraint 

\2 /10^^GpV\ 
m,, > 1.1 X 10-6 eV -- f{h/y/X, 0.8) = 7 x 10"^' eV , (28) 



h 
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which gives, together with P7|), a wide range of allowed m,^^ values. 

Finally, we want our RH neutrinos to decay before the Universe reheats via the decay of 
the Higgs, Tn-^ > = Hitj-ii). This implies a bound on the reheating temperature, Tj-h ~ 



o.i2Jr^Mp 13 , 



For Ml = 10^^ GeV, and chosing the mass of the light neutrino to be m,^^ = 10~^ eV, we find 
that the reheating temperature should satisfy < 2 x 10^° GeV, which is compatible with the 
bounds coming from gravitino production. 

Now we should estimate the effective temperature of the decay products of the RH neutrinos 
Ti, assuming that their energy density at the time of decay ti is converted into a thermal 
bath pNiiti) = ^g^:Tf. The neutrino energy density at ti can be calculated from its value at 
symmetry breaking as 

P.M={^)J^Y\*(h), (30) 



taking into account that, from tsB until ti, the RH neutrinos are non-relativistic and thus their 
energy density decays like matter, p_F ~ while the energy density in Higgs particles is 
dominated by the gradient and kinetic terms, with a radiation equation of state, ~ a~^. This 
gives the extra factor a(ti)/a(tsB) in ( pOD coming from the expansion of the Universe. Since 
at ti the energy density is dominated by the infiaton, we have p^iti) = 3M|r^y87r. For our 
parameters, we find 

Ti = 2 X 10^2 GeV 0^ X'/^'^ f^\h/VX, 0.8) (j^TK^) = 4 x 10^ GeV , (31) 

which does not constitute a problem with gravitino production. 

A final concern is that the lepton-number-violating processes among the particles in this 
thermal bath are always out of equilibrium, so that they do not wash out the lepton asymmetry 
produced, T^l=2 = ■:^Tf J^n^l^ < H{Ti). We have included here the possible channels of 
lepton- number violating interactions mediated by the heavy RH neutrinos N2 and N^. This 
bound gives 

Ti < 3.3 X 10^3 GeV ( , (32) 



which is easily satisfied, see Eq. (pTf), for the range of light neutrino masses consistent with 



neutrino oscillations [221. 



We may now ask how efficient is this mechanism for producing the required amount of 
baryons in the Universe, ti-q/ s = (4 — 9) x 10~^^, where s is the entropy density, and we use the 
fact that this ratio remains constant since reheating. 

The baryon asymmetry is produced via sphaleron transitions that violate (B+L) and convert 
a lepton asymmetry into a baryon asymmetry, = (28/79)nL ITHl' '^^^ entropy density at 
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reheating can be computed in terms of the energy density, s = (4/3)p<^/Trh, and can be estimated 
as 



s 



28 ul 
79 T 



rh 



28 STrh 

79 



4 



rh 



28 3nL r^jvi 
79 4n '^'^ 




A5/4 



2.2 X 10"^£— -/(VVA,0.8) 



h 



(33) 



where we have assumed that the lepton number density is directly related to the number density 
of RH neutrinos via nL = eni^^. For the values of the parameters we have chosen, the baryon 
asymmetry becomes u-q/s = 2.3 x 10~^£, which requires a relatively large leptonic asymmetry 
e = 3.4 X 10~^, that could be obtained in realistic neutrino models with large mixing angles. 
However, significantly smaller values of e can still be accommodated in this scenario since the 
values we chose for h and A were very conservative. If we keep the GUT vev at 10^^ GeV, we can 
still increase the production of RH neutrinos by taking the maximum values of the couplings 
allowed by the bounds (|2^ ) and (|^), h = 0.5 and A = 3 x 10~^, which gives e > 7 x 10~^ for 
agreement with the observed baryon asymmetry of the Universe. Moreover, the value of A can 
be chosen greater than the bound ( p4D if we accept that the Higgses thermalise (before decaying) 
soon after symmetry breaking, while keeping h small (PBD, so that the RH neutrinos remain out 
of equilibrium before decaying. 

We can also choose a larger vev for the GUT Higgs, e.g. v = 10^^ GeV, thus allowing a 
larger RH neutrino production. This vev value would be compatible with all the bounds, and 
permits a maximum value for A = 0.02 and h = 1.15, which gives ub/s = 1.12 x 10~^e. This 
corresponds to values of the leptonic asymmetry as small as e = 7 x 10^^, to be in agreement 
with the observed baryon asymmetry of the Universe. Smaller values of e would be difficult to 
accommodate in our scenario. 

We conclude that our scenario can generate a successful leptogenesis for a wide range of model 
parameters, although in order to go beyond our estimates, we would need more information on 
the structure of neutrino masses and mixings. 



5 Other relics, dark matter and cosmic rays 

This out-of-equilibrium production of particles during SSB may have other interesting cosmo- 
logical consequences. For instance, this mechanism may be responsible for the non-thermal 
production of superheavy dark matter (SDM) Consider the following scenario: a GUT 

scale symmetry breaking at the end of a period of hybrid inflation, where the SB field may 
generate via the mechanism described above, a large population of particles out of equilibrium, 
which may either decay into stable relics, or be stable themselves. Such a heavy relic will even- 
tually dominate the energy density of the Universe and constitute today the cold dark matter 
component, ficdm^^ = 0.14 ± 0.04, 95% c.l. where h = 0.72 ± 0.07 is the present Hubble 
rate in units of 100 km/s/Mpc. 
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Assuming that the rehc dark matter particles were produced soon after symmetry breaking 
and that the Universe reheated at a temperature T^h, we can estimate the fraction of energy in 
these non-relativistic particles today as 



^xh'' px T, 

rs^ 

Q^h'^ pSB To ' 



(34) 



where px = r py if the relic particles come from the decay, with branching ratio r, of species 
Y produced at symmetry breaking. Note that in the simplest case in which species X couple 
directly to the GUT Higgs, r = 1, but in general we expect r -C 1. Substituting the present 
temperature of the Universe, Tq = 2.73 K = 2.4 x 10~^^ GeV, we have 

px iq /lOi°GeV\ , , 

^ = 2 X 10-'' -— . (35) 

PSB V Jrh / 

Using expressions (0) and ( p^Sf ) , we see that in the limit of small couplings, px / Psb — 2 x 10~^ , 
where gx is the coupling to the Higgs, which gives it a mass Mx = gx^^'^ GeV. In order that 
the X-particle relics be the cold dark matter today, we require 



M^ = 10^GeV (i^^^j . (36) 



For a reheating temperature as high as Tj-h = 10^'^ GeV, the mass of the particle becomes 
Mx = 10'' GeV < Trh, and is thus non-relativistic at symmetry breaking. In this case, its energy 
density decays like radiation until its mass is of order the temperature of the Universe, so we 
should substitute Trh — > Mx in Eqs. ( p4D and (|35D. Therefore, Eq. ( p6D determines Mx = 10^ 
GeV, a very natural candidate for SDM. 

However, we must make sure that the SDM relics remain out of equilibrium since their 
production, otherwise their thermal population would be completely unacceptable today. For 
this we must ensure that their annihilation rate is always much smaller than the expansion rate, 
Tann = 'T-x (^Ih) CTann I'i' | "C H (Tj-i,) , where the annihilation crossection is cxann = 5'x/(167r^M|^). 
Using expressions ([T9|) and ( pO]) in the small coupling limit, nx = 4.8 x 10"'' A (^3^^ M|, we get, 
for A = 10-^, 

Fann = 3 X 10^ GeV g'^ « 200 GeV [j^^^) = H{T,^) . (37) 

It is clear that for the case above, with gx = 10~^, we are always safe. 

Alternatively, we may consider the production of X-particles from the decay, with very small 
branching ratio r ~ 10~^°, of very heavy Y-particles produced at SSB, with e.g. My = 10^^ 
GeV and Mx ~ 10^^ GeV. For this value of the coupling, gy = 10~^, the Y-particles remain 
always out of equilibrium, see (^^, while their energy density at symmetry breaking is given 
by py / Psb = 2 x 10^^, giving rise to a dark matter relic X which is compatible with present 
bounds. 
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Moreover, a very interesting consequence of this population of extremely weakly coupled 
heavy relics produced at GUT symmetry breaking is that they may constitute the origin of the 
ultra high energy cosmic rays |^5[. In the case the X-particles are metastable, with lifetime of 
the order of the age of the Universe, their mass Mx ~ 10^^ GeV = 10^° eV may be converted 
into extremely energetic particles that reach Earth today and whose nature could in the near 
future be studied by the Pierre Auger project, the high-resolution fly's eye, and the Japanese 
telescope array project |]25[|. 



6 Conclusions 

In this letter we have shown that the exponential growth of the Higgs vev towards its true 
vacuum at the end of a period of hybrid inflation, known as tachyonic preheating, can induce 
a significant production of particles, both bosons and fermions coupled to the Higgs. This new 
mechanism of particle production in the early Universe could be responsible for generating the 
observed baryon asymmetry via leptogenesis, as well as the present dark matter of the Universe. 

We have proposed a novel scenario for leptogenesis, in which a large population of out-of- 
equilibrium RH neutrinos is generated in the process of symmetry breaking after hybrid inflation, 
assuming that the symmetry breaking field is actually the Higgs field of a grand unified theory 
(GUT) with {B — L)-violating interactions. We have shown that this is indeed an attractive 
leptogenesis scenario, since it can explain the observed baryon asymmetry for a wide range of 
model parameters. 

In the case that stable relics are produced out-of-equilibrium from the decay products of the 
GUT Higgs, it is possible that they may constitute today the observed cold dark matter. These 
particles may have masses as large as 10^^ GeV, and could be responsible (if metastable) for the 
observed flux of ultra high energy cosmic rays. 
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Figure Captions 



Figure 1. The time evolution of the vacuum expectation value as compared with 

the approximate solution with mt* = 8. We have used a lattice of size N=128 and length 
L=1007r, which gives k^i^ = 2H = 0.02m and /Cmax = 2.22m. We also show the evolution of 
the number density, UBit), in units of 100 m~^, for bosonic particles coupled to the Higgs with 
(7 = 0.5. 

Figure 2. The spectrum of occupation numbers for both bosons and fermions in the (asymp- 
totic) true vacuum, using the lattice results for the Higgs vev, as compared with the analytical 
formulae (|T^ and (|T3|). The parameters chosen here are A = 0.01, g = 0.5 and h = 0.5. The 
tiny peaks at large momenta correspond to small resonances due to the strongly damped Higgs 
oscillations after SSB. 
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